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BACKGROUND OF THE INVENTION 

1. Related Application 

This application claims priority from co-pending provisional application Serial No. 
60/152,391 filed September 3, 1999. 

5 2. Field of the Invention 

This invention relates generally to the field of fluid dynamics. More particularly, the 
invention is a method for accurately and efficiently analyzing anisotropic turbulent flows of both 
gases and liquids. 

3. Background 

1 0 The study and engineering application of hydraulics or fluid mechanics dates to the dawn of 

civilization. The first irrigation canals were constructed before 5,000 BC. Theoretical analysis of 
fluid motions began in the modern scientific sense with Newton's Principia (1687) and his laws of 
motion. Bernoulli (1738) established a relationship between the pressure and velocity of an ideal 
fluid which became his famous theorem. d'Alembert (1744) deduced that the steady flow of an 

1 5 ideal fluid about a body produces no drag (d' Alembert's Paradox). Euler (1755) founded modern 
fluid mechanic analyses by deriving the differential equations of motion (continuity and momentum 
conservation) for an ideal (inviscid) isotropic fluid. Navier (1822) and Stokes (1845) 
independently derived the viscous terms necessary to extend Euler' s equations to isotropic, viscous 
fluids. The equation of state and the conservation of energy equation (First Law of 

20 Thermodynamics) written utilizing Fourier's (1822) law of thermal energy transfer completed the 
classical Navier Stokes equations governing isotropic fluid motions. 

Reynolds (1883) experimentally showed that there are two modes of fluid motion, laminar 
and turbulent motion, which occurs at large values of the dimensionless Reynolds number. 
Reynolds (1894), by considering a time average motion and a time dependent turbulent motion, 
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showed that turbulent motions introduce additional turbulent (Reynolds) stresses into the Navier 
Stokes equations, greatly increasing the number of flow properties required to describe turbulent 
flows. 

Despite over a century of intense aircraft and aerospace developments, wind tunnel and 
5 laboratory experiments, theoretical analyses and numerical computations based on the Navier 

Stokes equations, the determination of a general theory of turbulent motion remains one of the last 
great unsolved problems of classical physics. 

Maxwell (1858) founded modern kinetic theory when he introduced his new concept of the 
molecule velocity distribution function and deduced its equilibrium form. Maxwell (1867) derived 

1 0 the Maxwell transport equations and showed that the collision change integrals could be solved 
analytically without knowing the molecule velocity distribution function for Maxwell molecules. 
Boltzmann (1872) derived his integro-differential equation for Maxwell's molecular velocity 
distribution function and solved for his famous H-theorem, which proved that Maxwell's 
equilibrium form was correct. Chapman (1916), using Maxwell's transport equations, determined 

1 5 accurate general formulae for the gas transport coefficients. Enskog (1917) gave a general solution 
method for he Boltzmann equation, which showed that the Euler equations were the first and the 
Navier Stokes equations were the second approximate solutions to the Boltzmann equation and 
gave gas transport coefficients identical to Chapman's. Both the Chapman and Enskog analyses 
were for isotropic perfect gases, being based on perturbations of Maxwell's isotropic equilibrium 

20 solution. Their independent analyses cemented the belief in both the scientific and particularly the 
engineering professions that the Navier Stokes equations were generally and universally valid and 
that the correct way of analyzing both laminar and turbulent flows was through the Navier Stokes 
equations. This belief has been engraved in stone in both the scientific and engineering literature 
throughout the 20th century. 

25 Burnett (1935) derived the third approximate solution to the Boltzmann equations using 

Enskog' s solution method. Grad (1949) derived a thirteen moment method of solving 
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Boltzmann's equation which closely followed Enskog's methods. Both of these analyses have 
been proven to lack generality and have not been significant advances of the Navier Stokes theory. 
Bird (1963) introduced the Direct Simulation Monte Carlo (DSMC) numerical method for 
calculating ratified gas flows. Yen (1966) showed that the directional thermal energies 
5 (temperatures) were vastly different in shock waves and that the longitudinal temperature overshot 

(9 

its downstream value for shock Mach numbers greater than J— (1.34). Elliot and Baganoff 

(1974) showed that the Navier Stokes normal stress relationship and the Fourier energy flux 
component ratio were valid only in sound waves and were invalid at the shock end points. Elliot 

(1975) showed that the Navier Stokes normal stress relationship was invalid everywhere in shock 
1 0 waves and incompatible with the directional thermal kinetic energy moments of the Boltzmann 

equation. DSMC numerical shock calculations confirmed all of these analytical predictions and 
numerically illustrated the anisotropic not isotropic fluid features of many simple gas flows. 

Kliegel (1990), following a suggestion of Maxwell (1867) that near equilibrium flows were 
anisotropic and should be represented by an anisotropic, not isotropic Maxwellian, performed the 

1 5 first general anisotropic solution analysis of the Boltzmann equation. Kliegel (1990) showed that 
the Euler-Navier Stokes-Burnett equation sequence was not the correct approximate solution 
sequence for the Boltzmann equation. He showed that the correct gas dynamic equation set was an 
anisotropic fluid seven equation set for the density, three fluid velocity components and three 
directional thermal kinetic energies (temperatures), not the classic isotropic fluid five equation 

20 Euler-Navier Stokes set for the density, three fluid velocity components and total temperature. He 
also gave shear and directional energy flux relationship having the correct Mach number 
dependence for both sound waves and weak shocks. Kliegel' s anisotropic equations resolved 
d' Alembert's paradox giving a profile pressure drag associated with pushing an ideal (shearless, 
energy fluxless) fluid about a body. They also correctly predicted the directional thermal energy 

25 separations and overshoots and the thermal energy flux component changes occurring in shock 
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waves as predicted by Yen (1966) and numerical DSMC shock calculations. Bird (1994) 
summarizes recent DSMC calculational capabilities. Chen and Jaw (1998) presented a recent 
summary of classical isotropic fluid Navier Stokes based turbulent flow modelling. 

The present disclosure correctly consolidates all previous fluid dynamic and kinetic theory 
analyses and extends it to define a new method of analyzing and computing anisotropic turbulent 
flows. 

The prior art cited above may be found in the following references: 

1. 1687, Newton, L Philosophiae Naturalis Principia Mathematica, Oxford. 

2. 1738, Bernoulli, D., Hydrodynamica sive de Viribus etMotibus Fluidorum Commentarii, 

Berlin. 

3. 1744, d' Alembert, J., Traite de Vequilibre et du mouvement des fluides pour servir de suite au 

traite de dynamique? Paris. 

4. 1755, Euler, L., Principles genereaux du mouvement des fluides, Histoire de T Academie de 

Berlin. 

5. 1822, Navier, L.M.H., Memoires de V Academic des Sciences de Vlnstitut de France, V 6, p. 

389, Paris. 

6. 1822, Fourier, J., Theorie analytique de la chaleur, Paris. 

7. 1845, Stokes, G.G., On the Theory of the Internal Friction of Fluids in Motion, Transactions of 

the Cambridge Philosophical Society, V8, p. 287, Cambridge. 

8. 1858, Maxwell, J.C., Illustrations of the Dynamical Theory of Gases, Philosophical Magazine, 

V19, p.19 and V 20, p. 21. 

9. 1867, Maxwell, J.C., On the Dynamical Theory of Gases, Philosophical Transactions of the 

Royal Society, V 157, p. 49. 

10. 1872, Boltzmann, L., Weitere Studien iiber das Wdrmegleichgewicht unter Gasmolekulen y 

Styungsberichte, Akad. Wiss, Vienna, Part II, V 66, p. 275. 
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V 



11. 1883, Reynolds, 0., An Experimental Investigation of the Circumstances which Determine 

whether the Motion of Water Shall Be Direct or Sinuous, and the Law of Resistance in 
Parallel Channels, Philosophical Transactions of the Royal Society, V174, p. 935. 

12. 1895, Reynolds, O., On the Dynamical Theory of Incompressible Viscous Fluids and the 

5 Determination of the Criterion, Philosophical Transactions of the Royal Society, V186A, p. 

123. 

13. 1904, Prandth, L., Uber FlUssigkeitsbewegung bei sehr kleiner Reibung, Proceedings of the 

Third International Mathematical Conference, Heidelberg, Leipzig. 

14. 1916, Chapman, S», The Kinetic Theory of Simple and Composite Monotonic Gases: Viscosity, 
1 0 Thermal Conduction, and Diffusion, Proceedings of the Royal Society London, A, V 93, 

p.L 

15. 1917, Enskog, D., Kinetische Theory der Vorgange in massig verdunnten Gasen, Almquist & 

Wilsells, Uppsala, Sweden. 

16. 1935 Burnett, D., The Distribution of Molecular Velocities and the Mean Motion in a Non- 
1 5 Uniform Gas, Proceedings London Mathematical Society, V 40,p. 382. 

17. 1949, Grad, H., On the Kinetic Theory of Rarefied Gases, Convention of Pure and Applied 

Mathematics, V 2, p. 331. 

18. 1963, Bird, G.A., Approach to Translational Equilibrium in a Rigid Sphere Gas, Physics of 

Fluids, V 6, p. 1518. 

20 19. 1966, Yen, S-M, Temperature Overshoot in Shock Waves, Physics of Fluids, V9, p. 1417. 

20. 1970, Chapman, S., Cowling, T.G., Mathematical Theory of Non-Uniform Gases, Third 

Edition, Cambridge University Press. 

21. 1974, Elliot, J.P., Baganoff, D., Solution of the Boltzmann Equation at the Upstream and 

Downstream Singular Points in a Shock Wave, Journal of Fluid Mechanics, V 65, p. 603. 
25 22. 1975, Elliot, J.P., On the Validity of the Navier-Stokes Relation in a Shock Wave, Canadian 
Journal of Physics, V 53, p. 583. 
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23. 1990, Kliegel, J.R., Maxwell Boltzmann Gas Dynamics, Proceedings of the 17 th International 

Symposium on Rarefied Gas Dynamics, Aachen 1990, Edited by Alfred E. Beylich, VCH, 
New York. 

24. 1994, Bird, G.A., Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Oxford 

University Press, New York. 

25. 1998, Chen, C-J. Jaw, S-Y, Fundamentals of Turbulence Modeling, Taylor & Francis, 

Washington. 
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SUMMARY OF THE INVENTION 

This disclosure concerns a simple method of deriving the correct kinetic anisotropic fluid 
turbulent flow equations from the Boltzmann equation and the correct thermal and turbulent 
moment closures to yield a new, closed set of fluid dynamic equations describing anisotropic 
5 turbulent motions, solving the century old quest for the solution to this classical problem. 

Current turbulent flow analyses are based on the isotropic Navier Stokes equations, turbulent 
fluctuation moments of these equations and modeling of the unknown turbulent moments, identifying 
introduced unknown modeling coefficients by experiment. This results in a variety of semi-empirical 
turbulent flow analyses based on different Navier Stokes moment relationships satisfying mass, 
1 0 momentum and total energy conservation for an isotropic fluid, calibrated to approximately match a 
limited number of experiments. The limitations of current Navier Stokes based turbulent analyses are 
well known and appear to be irresolvable after a century of intense analytical and experimental effort. 

The present invention provides a general, closed, anisotropic kinetic turbulence theory for gases 
and liquids based on new solutions of the Maxwell moment equations of the Boltzmann equations. 

1 5 These solutions provide a closed initial equation set for the four time average fluid mechanic variables, 
the sixteen time average thermal motion correlation and the sixteen time average turbulent motion 
correlations listed in Table L Higher order closure sets involving higher order (fourth and fifth, etc.) 
thermal and turbulent moment sets are easily generated allowing solutions of turbulent flow problems 
to higher and higher accuracy if desired or needed for special problems. These higher order sets will 

20 not be further discussed since they add only unnecessary detail to the disclosure and reveal little that is 
new. 

An algebraic differential equation solution method described in my co-pending application 
entitled "Method for Algebraically Solving Differential Equations, including Stiff Equations, to High 
Accuracy" may be used to reduce the thirty six governing differential equations to the solution of as 
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few as six differential equations (continuity, momentum, total energy and turbulent energy) plus thirty 
algebraic relationships. The actual equation set solved depends on the problem and accuracy desired. 

The presented analytic anisotropic turbulence theory is mathematically correct, involving no 
semi-empirical modeling. The even order higher order moment closures involved are general and 
correct to known accuracy , The analysis can be used to solve the infinite set of Maxwell moment 
equations to any level desired to verify the accuracy of the initial closure set. 

The anisotropic kinetic turbulence theory presented allows the calculation of turbulent motions 
with the same accuracy and computer resources as current laminar flow calculations. 

The new teachings of the present invention may be summarized as follows: 

1. A closed turbulence theory may be simply derived directly from the Boltzmann 
equation without turbulent correlation modelling or experimentally determined 
modelling coefficients. 

2. Turbulent flows are totally inviscid flows, unchanged by molecular collisions. 

3. Turbulence is universal, having the same governing equations and the same physical 
structure in both gases and liquids because turbulence is an inviscid flow 
phenomena described entirely by inviscid convective flow equations, independent of 
collisions and molecular interaction laws. 

4. The fact that time average turbulent flow correlation between turbulent flow 
fluctuations and density or thermal velocity correlations are zero because there is no 
physical (collisional) connection between these fluctuations, since both the turbulent 
flow Auctions and density are conserved in collisions. 

5. A method of analyzing the Maxwell moment equations to obtain a closed set of 
equations which describe time average turbulent flows to any desired accuracy level. 
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All higher moment equation above the closure level can be solved in even/odd sets of 
two, thus yielding solutions to all Maxwell moment equations with error of Burnett 
order or less. 

Direct predictions of turbulent transition locations, turbulent energy growth rates, 
shears and other thermal and turbulent flow quantities of physical importance can 
now be derived or calculated from solution of the presented equations with the same 
accuracy and computer resources as current laminar flow analyses or calculations. 

Universal even order thermal and turbulent moment closure relationships exist and 
are explicitly given, allowing the termination and solution of the Maxwell moment 
equations at any desired odd moment equation set, and solution of the Boltzmann 
equation to any desired approximation. 
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DETAILED DESCRIPTION OF THE INVENTION 



In the following description, for purposes of explanation and not limitation, specific details 
are set forth in order to provide a thorough understanding of the present invention. However, it will 
be apparent to one skilled in the art that the present invention may be practiced in other 
5 embodiments that depart from these specific details. In other instances, detailed descriptions of 
well-known methods and devices are omitted so as to not obscure the description of the present 
invention with unnecessary detail. 

The Boltzmann equation may be written in fluid dynamic collision rate form: [See 
Appendix, Equation set 1]; 

1 0 where ^ = [pF(v)j is the Lagrangian convective motion derivative, p is the fluid density, t> is the 
molecular collision rate, F(v) is the molecular velocity (v) distribution function and F(v)* is the 
molecular velocity distribution function collided into F(v). There is no approximation in writing the 
Boltzmann collision integral in this format, since the required thermal velocity moments of F(v)* 
are correctly calculated from the full Boltzmann collision integral (Kliegel (1990)). The molecular 

1 5 collision rate is 

V = p K 

for Maxwell molecules, where K is a molecular collision constant. 

The illustrative derivation will be done for Maxwell molecules for simplicity, since the 
collision change integrals are known exactly for Maxwell molecules. The same derivation can be 
performed for any molecular model, it is just more algebraically complex and adds nothing to the 
20 teachings. 

Both the original Enskog (1917) solution as a perturbation about the equilibrium 
Maxwellian and the Kliegel (1990) solution as a perturbation about a non-equilibrium anisotropic 
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Maxwellian are Navier Stokes order accurate ( —r Burnett order error terms) and do not depend on 

v 

density gradients. Thus, density gradient corrections to the molecular velocity distribution function 
and its moments belong to the Burnett order of corrections. This important fact, although present in 
the result of both analyses, is explicitly stated by neither author nor utilized to further extend or 
5 simplify the analyses. 

Multiplying the Boltzmann equation by the molecular velocity product vjv™ V3, where /, m, 
n are integers, and integrating over all possible velocities (indicated by the bar) yields the Maxwell 
molecular velocity moment equation: [See Appendix, Equation Set 2]. Derivation of the 
anisotropic turbulent flow equations from this equation will be discussed later. 

1 0 Consider the molecular velocity (v) to be composed of a mass mean fluid flow velocity (u) 

and a thermal motion velocity (c), then the Cartesian components of these velocities are related by 

Multiplying the Boltzmann equation by the thermal velocity product cjc™ C3 , where /, m and n are 
integers, and integrating over all possible thermal velocities (indicated by the bar) yields the 
1 5 Maxwell thermal velocity moment equation: [See Appendix, Equation Set 3], 

For odd moment (/ + m + n = 2N + 1, N integer) equations, the coefficients of the density 
gradient terms are even moments and products of lower even moments, which have equilibrium 
values. Thus these even moment coefficients must be zero to Burnett order for the density gradient 
terms to be of Burnett order. This requirement yields the general even order (fourth, sixth and 
20 higher) density gradient independent thermal moment closure relationships: [See Appendix, 
Equation Set 4]. 

The fourth order thermal moment closure relationships can also obtained from Kliegel 
(1990) by simple calculation but he did not state these results nor utilize them in any manner. 
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Substituting these relationships into Equation Set 3, one obtains the closed odd order (third, fifth 
and higher) density gradient independent thermal moment equation. [See Appendix, Equation Set 
5]. 



The classical Maxwell thermal velocity twenty moment equation set is: [See Appendix, 
5 Equation Set 6]. 

After using the fourth order moment density gradient independent closure relationships: 
[See Appendix, Equation Set 7], the closed directional thermal energy flux and structure equations 
become: [See Appendix, liquation Set 8] where the fourth order thermal moments have been 

* 

replaced by products of second order moments through the fourth order density gradient 
1 0 independent closure relationships. 

The above set of equations are the correct anisotropic gas dynamic flow equations with 
neglected terms of Burnett order in the convective terms of the closed thermal energy flux and 
structure equations. All other equations in the set are exact. These equations resolve the known 
deficiencies of the Euler-Navier Stokes-Burnett equation sets, d'Alembert's paradox, the normal 

1 5 stress incompatibility, the incorrect Mach number dependency of the Navier Stokes normal stress 
and Fourier energy flux components, and the incorrect coordinate rotation variance of the Burnett 
equations. These equations are valid for both sound waves and weak shocks, and thus will 
accurately predict all flows except in those regions of high change rates in distances comparable to 
the mean free path length, such as occur in strong shock waves and at boundaries in very rarified 

20 flows. 

The twenty moment equations may be further simplified by performing an Enskog-type 
ordering analysis on all terms to obtain the reduced equation set: [See Appendix, Equation Set 9] 
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One can extend the above analysis to time average turbulent flows by considering the mass 
mean fluid flow velocity (u) to be composed of a time average velocity (u) , and a turbulent motion 
velocity (u')> like Reynolds (1894). The Cartesian components of these velocities are related by: 

u x = u x + u[ 

5 where now 

The turbulent motion (u') is inviscid, unaffected by collisions since the mean flow velocity 
(u) is conserved in collisions and only the thermal (c) velocities are changed in collisions. This 
important fact has not been previously recognized, explicitly stated or utilized in the 

1 0 scientific/engineering literature. Since both the density and turbulent motion are preserved in 

collisions, there is no collisional (physical) coupling between turbulent motion fluctuations, density 
fluctuations and thermal velocity correlation fluctuations. Thus, there is no time average correlation 
between these fluctuations. The averaging time required for the time average correlations between 
turbulent velocity fluctuations and molecular (number density and thermal motion correlation) 

1 5 fluctuations to be essentially zero is quite short, since the molecular fluctuations occur on a 

collisional time scale of nanoseconds and turbulent flow fluctuations occur on much larger fluid 
dynamic time scales (typically measured in milliseconds or greater). Thus, the time dependent time 
average anisotropic turbulent flow equations presented can be used to compute turbulent flows 
changing on fluid dynamic time scales for the vast majority of engineering flows. The time change 

20 adequacy of the time average anisotropic turbulent flow equations needs to be verified for extremely 
rapid flow changes such as those associated with shocks, flames and deteriorations. The only 
turbulent flow variables necessary to analyze anisotropic turbulent flows in the vast majority of 
physical cases are those given in Table I below, all other time average correlations being zero. 
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Table I 
Turbulent Flow Variables 



Fluid Dynamic 


Thermal Moments 


Turbulent Moments 


P 


(i) 




(3) 


< 


(3) 




(3) 




(3) 


uX 


(3) 








(3) 




(3) 






C 1 C 2 


(6) 




(6) 






CjC 2 C 3 


(1) 


u;u^u; 


(1) 


Four Fluid Dynamic 
Variables 


Sixteen Thermal Moment 
Variables 


Sixteen Turbulent Moment 
Variables 



The time average turbulent flow variables are the flow density (p) , where the bar indicates 
5 time averaging, the three mass mean velocity components (u x ), the three anisotropic thermal kinetic 
energies ( c \ , where the two bars indicate both time averaging and averaging over all thermal 
velocities), the three thermal shear components (c^), the nine thermal kinetic energy fluxes (Cjcf 
and c x c 2 2 ), the thermal structure correlation (c 1 c 2 c 3 ), the three anisotropic turbulent kinetic energies 
( u[ 2 ), the three turbulent shear components ( u^ ), the nine turbulent kinetic energy fluxes ( u^ 2 
1 0 and u^ 2 ) and the turbulent structure correlation ( uju^ ). 

The instantaneous Maxwell molecular velocity (v) and the thermal velocity (c) moment 
equations (Equation Sets 2 and 3) are both valid at every instant of time and the instantaneous 
Maxwell thermal velocity moment equations can be derived from the instantaneous Maxwell 
molecular velocity moment equations. The time average Maxwell molecular velocity (v) and time 
1 5 average thermal velocity (c) moment equations are also both valid simultaneously. The time 

average thermal velocity equations yield equations governing the behavior of the anisotropic time 
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average thermal directional kinetic energies, shears, directional kinetic energy fluxes and structure 
correlations: [See Appendix, Equation Set 10]. 

The time average molecular velocity equations yield equations governing mass and 
momentum conservation and the behavior of both the time average thermal and turbulent 
5 directional energies, shears, energy fluxes and structure correlations: [See Appendix, Equations 
Set 11]. 

After eliminating the collision change terms with the time average thermal Maxwell 
moment equations and other simplifications, one obtains the equations governing the behavior of 
the anisotropic time average turbulent directional kinetic energies, shears, directional energy fluxes 
1 0 and structure correlations: [See Appendix, Equation Set 12]. 

Both sets of third order directional kinetic energy fluxes and structure correlation equations 
are closed using the fourth order density gradient independent time average thermal and turbulent 
moment closure relationships: [See Appendix, Equation Sets 13 and 14]. 

One thus obtains the complete closed Maxwell twenty moment equation set governing 
1 5 anisotropic time average turbulent flows: [See Appendix, Equation Set 15]. 

By the same process, corresponding higher order time average turbulent moment equation 
sets may be derived using higher order closure relationships. The general density gradient 
independent time average thermal and turbulent even moment closure relationships are: [See 
Appendix, Equation Sets 16 and 17]. 

20 These closures allow the infinite set of time average Maxwell moment equations governing 

turbulent flows to be closed at any desired odd order level starting at the third order. If one closes 
at the third order (say), then one can solve all higher order moment sets (fourth and fifth, six and 
seventh, etc.) in order, where the errors are still of Burnett order in the convective terms of the 
highest order moment equations (fifth, seventh, etc.) considered in the set. The method thus allows 
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solution of all the Maxwell moment equations, both those sets below and above the chosen closure 
level. 



The above results may be simply extended in many ways. For example, the general even 
order density gradient independent moment closure relationships (Equation Set 4) are valid in 
5 equilibrium. Thus, [See Appendix, Equation Set 18] and the even order collision integrals 



It will be recognized that the above-described invention may be embodied in other specific 
forms without departing from the spirit or essential characteristics of the disclosure. Specifically, 
the anisotropic turbulent flow equations may be reformatted and/or approximated and computed in 

1 0 many different ways to solve specific problems. Obvious first steps would be to add gravitational 
forces or to sum the directional turbulent energy equations to obtain a total turbulent energy 
equation or to reduce the complexity of the equation set through a Prandtl (1904) boundary layer 
analysis. The analysis presented (valid for monatomic perfect gases) can be expanded to diatomic 
(and polyatomic) gases, gas mixtures, reacting gas mixtures, plasmas, etc. through known methods, 

1 5 Chapman and Cowling (1970). The analysis can be extended to isotropic fluids by setting the three 
directional thermal kinetic energies cf or cf equal, and solving only the summed total thermal 
energy equation, rather than the three directional thermal kinetic equations. Since this reduces the 
Equation Set 9 to the Navier Stokes equations, Kliegel (1990), the analysis can be extended to 
liquids by setting both the three directional thermal kinetic energies equal and the density constant. 

20 The closed anisotropic turbulent flow equation set can be utilized with existing Reynolds averaged 
Navier Stokes equation solvers to solve incompressible isotropic fluid problems. Thus, it is 
understood that the invention is not to be limited by the foregoing illustrative details, but rather is to 
be defined by the appended claims. 
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c i + c 2 c 3 ~ c i + c 2 c 3 > etc - ma Y be calculated for all molecular models. 



While all equation sets may be solved by straightforward computation, they may be more 
efficiently solved using the method described in my co-pending application entitled "Method for 
Algebraically Solving Differential Equations, Including Stiff Equations, to High Accuracy". 
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APPENDIX 



Equation Set 1 
(1-1) 

£[PF(V)] = ^pF^.v.Af.p^.^^pp^^^^P^ 
= pt>[F(v)*-F(v)] 
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Equation Set 2 
(2-1) 



fHv 2 m v 3 "] + A[ pv ; 



p v v;v-v 3 " 



■ vVV 

v l v 2 v 3 



+ ■ 



-[pvivrv-]- 



d r 



where v(v™ v£ is the instantaneous molecular velocity moment vjv^ V3 averaged over F(v) and 

* ^ 

v i v 2 v 3 is tlie instantaneous scattered molecular velocity moment vjv™ V3 averaged over F(v) . 
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Equation Set 3 
(3-1) 







<2k. 


12 3 



/Cj c 2 c 3 



+ i[ c ' Ci ) + ^[ c ' c3 ] 



- /2C X C 2 C S 



y A i v/v 2 ^3 



+ / 



4- m 



12 3 1 2 



3 

c M c m-I n UKX 2 , / ^ u 2 , ~ 0 T m-/ n¥l <^2 

c l C 2 *"3 ^ L 1 C 2 C 3 -x ^ U 1 C 2 C 3 "^T" 

esq chi 2 chi 3 



+ n 



eta. 



3 



+ 



c[ +/ c^»/(cf) c^cj-z^cj cic^-^^Cg] cf<c 
cfc^c'-zfc^) cpc^-^[c|] C^Cj"^^] 

/f3 c^c 3 *-«[Sa] cfcrcj-^p] cfc-cr 



3 



+ icc-cr 7 



* n 

C 1 C 2 C 3 _C 1 C 2 C 3 



1 8p 

1 3p 

p3x 2 

1 3p 

- p 3x 3 



where cjc^cj is the instantaneous thermal velocity moment cjc^ averaged over F(v) , and 
* 

c[c 2 c 3 is the instantaneous scattered thermal velocity moment cJc^Cj averaged over F(v)*. The 
latter may be calculated from the full Boltzmann collision change integrals, Kliegel (1990). 
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Equation Set 4 
(4-1); (4-2); (4-3) 



cf 7 ^ = /[cf] c{- 7 c*c^+i»[c^] c|cr^+n[^] ^^f 7 
clefts = /[^] c^c^ + mp] c{cr^"+npS] 



neglected terms being of Burnett order, where /+m+n = 2N +l,an odd number. 
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Equation Set 5 
(5-1) 



d_ 

dt 



cjc-c;] + Ui £-[c[c^] + u 2 ~[o[c m 2 c" 3 ] + u 3 ^-[cjc^] 



d_ 



+ / 



+ m 



n 



= v 



W ^2 ^3 -i ^ C 1 C 2 L 3 ^^lW C 3 n 

ox 1 dx n chi 3 



<2k 2 



„n-7 ^3 , „l„m+l„n-l ^3 



<2u. 



: 1 C 2 C 3 C I C 2 Zc 3 j 



where /+m+n = 2N +1, an odd number and the neglected terms are of Burnett order in the 
convective terms (left hand side). 
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Equation Set 6 
(Twenty Moment Equations) 



Continuity (6-1) 



Momentum (6-2) 

f [ p + 8^[ P t U ' + °' 1] + 8^[ p l UlU2 + + ^"t P I UlU3 + ^JH 
Directional Thermal Energy (6-3) 

lP] +u 'S:R] +u ^[ 



_3_r- 

dx. 



R] 



+ 2 



— du, — du, — 3u, 

C '- +C ' C ^ +C ' C 'ST 



3xj 



Thermal Shear (6-4) 

at L ^ J ' " x ax 



^[c t c 2 ] + u x ^-[^cj + u 2 ^[c^] + u 3 ~-[c lC2 ] 



dx. 



t 3u, — 3u, du, 

+ C 1 C 2 T^ + C 2 ^ + C 2 C 3 



3xj 



dx 0 



dx? 



T 3u 9 diu du,, 



+ <r 



i 

+ — 
p 



■H-CX 



■ + c,c 



3xj 1 2 dx 2 1 3 dx^ 



— 1) CjC 2 
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Equation Set 6 

Directional Thermal Energy Fluxes (6-5) 
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(6-6) 



Equation Set 6 



3 r — 2I d 





d 


Clc 2. 


+ u 2 — 


dx, 



c f c;i + u 3 — ic^ 



dx. 



— 2 — 2 du t — 2 ^Ui 
ox, dx 0 dx. 



+ 2 



- c 0 





du 2 


C 2 C! 


dxj 


" a 




dX( 





+ C 1 C 2 ^ 2 +C 1 C 2 C 3 ^ 

dx 2 dx 3 



+ ^-[c 1 c 2 ] + ^-[c 1 c 3 ] 



3x 0 



2 CjC 2 



d r— n | d 

dx^ 1 2j 3x 2 



+ 



9xj 



dx, 



+ [c?c*-c? [c2]-2c,c 2 [c 1 c 2 ]] 



1 dp 
p9x, 



+ |c,c 2 c 2 -3 0,02 



— r-Iiil dp 



p dx 2 

C1C3C2 - [4] ~ 2 Cfi~ 2 ff" 



= \ [Jc 2 c 2 -3c 2 cf] 
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Thermal Structure (6-7) 



Equation Set 6 



-[ Cl c 2 c 3 ] + u t ^-[cic 2 c 3 ] + u 2 _[ Cl c 2 c 3 ] + u 3 3^-[c,c 2 c 3 ] 



du, — 2 3Ui — 2 



+ C 3 C 1 



2 



+ C 1 C 2 C 3 



3u 0 



+ C 0 Q 



dx r 

2 5U 3 ~T2 du 3 



3x, 



• + c T c: 



1^3 



2 V 1 



c 2 c 3 



3x, 



+ c,c 



l w 2 



dx, 



■ + c 1 c 2 c 3 



2 d}h 
3u 3 



_3_r- 
3x, 



Kl + ^[ CiC2 l + ad CiC3 l 



3x, 



- C 1 C 3 



Cl c 2 



3xj 



9x, 



3x, 



[c 2 c 3 c 2 ] + ^-[c,c 3 c^] + A[ CiCaC 2 j 



3x, 



+ [c 2 c 3 c 2 - c^" [of] - 2 
+ [c,c 3 c 2 - [c 2 ] - 2 
+ [c^c 2 - [c 3 2 ] - 2 [c^]] 1 df> 



p 3x 3 



- V CjC 2 C 3 



where c 2 =cf + c 2 + c 2 and the collision change terms (right hand sides) have been calculated 
using the method of Kliegel (1990). 
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Equation Set 7 
(7-1; 7-2; 7-3; 7-4) 



2 2 q 2 L 21 
C 1 C 1 = 3C 1 [ C lJ 

Cj^cf = 3 [cf ] 

c i c 2 c 3 = [c*l + 2 [c^] 
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Equation Set 8 
Closed Directional Thermal Energy Fluxes (8-1; 8-2) 



3 r — 51 d 
+ 3 



ClC '] +U2 aB ClCl2 l +U3 ^:[ ClC '] 



a_ 

dx. 



2 ^Hl 

C[Cl dx, + ° 2Cl 3x 



2 3Uj ^-^2 ^1 



+ C 3 C 1 V- 
2 5x 3j 



+ 3 



3x, 



[c 2 ] + Cl c 



= |[ c i c2 -3c lC f] 



3x, 



3 r-7i 



+ C l C 3-^— l C i 



3x, 



a r 2 

— C-C 

at L 



3 r — 2I ^ r — "i 

cxJ + m 



2 

2 



2 2 3u x 

ox j dx 2 dx 3 



+ 2 



3u 



3u 9 



C 2 C 1 2 +C 1 C 2 

dXj dx 2 



+ CjC 2 C 3 



3u 2 
3x. 



+ Of ^-[C^] + 0,0, ^-[cj] + 0,03 



9x 1 



9x 0 



9x. 



+ 2 



3 r 1 — d r -] 3 

C 1 C 2 — [ C 1 C 2 J + C 2 ^"l C l C 2 J + C 2 C 3 



9xj 



3x, 



[c,c 2 ] 
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Equation Set 8 

Closed Thermal Structure (8-3) 

3u 1 — J 3ui — 2 3u, 

+ W3 T^ + ^c, t-^ + c^ — L 
dXj dx 2 dx 3 

— j 3u 2 d\x 2 — j 5u 2 



+ c 3Cl 3xi +W3 9x2 +¥ 3aXj 

2 9u 3 2 ^ U 3 ^U 3 

+ C 2 C 1 ^ '" C 1 C 2 ^ l " C l C 2 C 3 ^ 

aXj dx 2 dx 3 

+ C l C 3^-[ C l C 2] + C^^-[cIc7] + C3^-[clc^] 



the neglected terms are of Burnett order in the convective terms (left hand sides). 
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Equation Set 9 



Continuity (9-1) 



Tw + 4tIp U J + ^-[p u 2 ] + ^-[p u 3] = 0 



3t 3xj 



3x, 



Momentum (9-2) 



3u, du, 3ui du, 
— 1 + u 1 — 1 + u 2 1— 1 + u 3 — — 



at 



= 0 



3xj 3x 



3x, 



1 


" a 




p 


3x 1 


p cf 



Directional Thermal Energy (9-3) 



at 

+ 2 
1 



auj 



duj 



3uj 



P 



c, — h CjC 2 3-^ + c,c 3 
ox, dx 2 dx 3 

ax 2 '- r ' ~' L ' li+ Bx 3 



_a_ 

ax, i 



[pc^J + ^p c 2 cf] + ^-[pc 3 cf] 

3C -c t j 

Thermal Shear (9-4) 



<>iC 2 = — 



_1 



~r 3u< — 3u 



3x 9 



3x x 



Directional Thermal Energy Fluxes (9-5) 

2v[ dx^ 1 ! ax,L 2 J a^L 3 J 



2 Ci 
c lCl = l - 



* ^ 2 ^2 
C 2 C 1 ~ 

21 6v 



Equation Set 9 



9 r^i . 9 rr 



3x, 



3x 2 L 
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Thermal Structure (9-6) 
C!C 2 c 3 = 0 

neglected terms being of Burnett order in the thermal shear, directional energy fluxes and structure. 
The shear stress and thermal energy flux relationships are identical to the results Kliegel (1990) 
obtained by an Enskog-type solution of the Boltzmann equation. 
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Equation Set 10 
Time Average Directional Thermal Energy (10-1) 



d_ 

at 

+ 2 



1 
P 



+ u 



1 dx, L 1 



- a 



9x 0 



1 3x x 



l v 2 



3Uj 



-hex 



3u, 



3xo 



l v 3 



2 3x 3 



3xj 1 1 J 3x 2 L 



3 r- — 

— p c^ 



+ ■ 



3x, 



Time Average Thermal Shear (10-2) 

a r==1 _ a , _ a 



|[c x c 2 ] + u, ^-[ Cl c 2 ] + u 2 ^-[ Cl c a ] + u 3 ^-[ Cl c 2 ] 



= 3u, -j 3u, = 
dx t dx 



2 

air 



2 ^ X 3 



= 3u 9 = 3u, 
+ c, "T^ + CiC, 3-^- + c x c 3 — 
ox, dx, 



3x, 



1 

+ = 
P 



" a 




c 


' r- = =i 






axj 


pc 2 c* 


+ a* 


[p c i c 2 j 


+ ax~t p ClC2 ° 3 - 
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Equation Set 10 



Time Average Directional Thermal Energy Fluxes (10-3) 



d 

3tL 



c i c i 



-JL 

3x, 



JjCf j + u 



- a 



3x, 



a r 



+U3 ax7l ClC v 



+ 3 



+ 3 



Cl c, 



3uj 
dxj 



+ C 2 C 1 



2 <K , 73 3^ 



• + c,c 



3x 2 3x 3 



1 <3xj 



+ Cl c 2 



= a 
+ c i c 3 37 



2 



2 o 2 
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Equation Set 10 

Time Average Directional Thermal Energy Fluxes (10-4) 



d_ 
3t 



+ u~— 
1 dx 1 



— d 

+ U, 



dx 2 i 



Cl c 2 2 



+ u. 



== 9u, == 3u, =f 3u, 



3xj 



3x, 



3x, 



+ 2 



=f 9u, =f 3u, ===== 3u 

C 2 C, rr-^- + C,C 2 — ■ i + C 1 C 2 C 3 



9xj 



3x 2 9x 3 



+ 2 



+c: 



<2x, 



c,c 2 ] + c 2 c 3 ^-[ciC 



^N] +CiC2 £Nl +CiC3 ^:N] 



d 



1 2 _ o 2 



„ 2 V,, 
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Equation Set 10 



Time Average Thermal Structure (10-5) 




FILE NO. 004210.P001 -35- GWH\DLM 



Equation Set 1 1 



Continuity (11-1) 



Momentum (11-2) 

|[p^]^[p[uT 2 + < + ?]] 

+ ^-[p[u7u7 + <ul + c^]] 



+ ^-[p[uiU 3 +<u;+C lC 3 

= 0 



Turbulent Directional Energy (11-3) 



dt 



[p [u, 2 + u[ 2 + c 2 ]]+^-[p [u/ +3 u,u[ 2 +« 2 + 3u lC 2 +c lC 2 ]J 

-}] 
-}] 



d r- r — 2 . — —pi 



9 r- r — 2 — — 7T 



P | U 3 U / +U 3 Uj Z -h2u 1 UjU3+U3u[ 2 + U 3 Cj + 2 U 1 C 1 C 3 + C3C 



pv 



3 1 
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Equation Set 1 1 



Turbulent Shear (11-4) 



3t 



u 2 Uj 2 +u 2 u[ 2 +u 2 u( z +2u x u[u 2 +u 2 Cj +2u 1 c 1 c : 



]] 



d r- r — 2 . — — 72 



+ ^"[^[ U 1 U 2 +U 1 U; Z +2a 2 U;U2+<U; +2U 2 C!C 2 +11^2+0; 

= - pV CjC 2 

Turbulent Directional Energy Fluxes (11-5) 



J""[p [ u i 3 + 3 u i u i 2 +U X 2 + 3u i c ? +c i c ?]J 

+ irfi P +6 ^ +4 <v ]] 

+ ^- j^p ^6 u7 2 cf + 4 u7 c t Ci + 6 u[ 2 cf + c 2 c 2 j j 

+-^-[p [u7 3 u7+3u7 2 u]ul + 3u7u 2 "u; 2 +u7u 1 '< 2 +3u7u 2 u 1 ' 2 + u 1 'u 2 u; 2 ]] 



a r-r„-2 — 



|p| Su/c^ +3 UiU 2 cf+u 2 c 1 c 2 + 3u!C 2 Cj 2 +3u[ z CjC 2 +3u 

3x a 



^cf + c^c 2 ]] 



[p [u t 3 u 3 +3u! 2 u'u3 +3u 1 u 3 u 1 ' 2 +u 3 <uf +3u x u^ + u^uf ]] 



+-^-| p| 3u 1 'c 1 c 3 + 3 UiU 3 Ci +u 3 c 1 c 1 z +3u 1 c 3 cj I +3u' z CjC 3 +3u^u 3 c 2 +c t c 3 c z 

OX -3 









pv 


3uj 


3 1 



1 

+- 

2 



2 o 2 
CjC -3 CjCj 
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Turbulent Directional Energy Fluxes (11-6) 



3 f- T 2 . 77 



+ ^—[p [ ^i 2 ^ 2 + ^ 2u 2 2 + ^2 2u i 2 +4u^u^u[u2 + 2u^u 2 2 + 2u^u 2 u[ 2 + u[ 2 u 2 2 jj 



a r-i 



+ ^~[p[ u i 2c 2 +u 2 2c2+ 4u 1 u 2 c 1 c 2 +2u 1 c 1 C2+u[ 2 C2+2u 2 c 2 c 1 z + U2 Z c 1 z +4u[u2C 1 c 2 +C1C2 

r- r — 

dx 2 

+ "^~[^[ 3u 2 2 c 1 c 2 + 3 UjU 2 c 2 +Uj c 2 c 2 + 3u 2 CjC 2 +3u 2 2 c 1 c 2 +3u^u 2 c 2 -t-CjC^jj 

+-^-[p [ u 1 u 3 u 2 2 +u 1 u 3 u 2 2 + u 2 2 UjU 3 +2 UjU 2 u 2 U3 +2 u 2 u 3 u[u 2 +2u 2 UjU 2 u 3 + UjU 3 u 2 2 +u 3 u 1 'u 2 2 + u 1 'u 3 u 2 2 
dx. 1 1 



.'2^2 



2/1 



[p [ u i u 2 3 +3u 2 2 UjU 2 +3u 1 u 2 u 2 2 +u 1 u 2 u 2 2 +3u 2 u[u 2 2 + u[u 2 u 2 2 Jj 



+- 



3x 3 L 



u 2 2 c 1 c 3 + u 1 u 3 c 2 +2u 1 u 2 c 2 c 3 +2U2U3CJC2 +UjC 3 c 2 +2u 2 CjC 2 c 3 + u 3 c 1 c 2 jj 
u 2 2 CjC 3 +2UJU2C2C3 + u[u 3 c 2 +2u 2 u 3 c 1 c 2 +CjC 3 C 2 Jj 



.2 — 



= pv 



u. 



-2u 2 c 1 c 2 + 



1 



CjC 2 -9 c x c. 
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Turbulent Structure (11-7) 



dt 



P [u 1 U 2 U 3 +U 3 C 1 C 2 + U 2 C 1 C 3 +U 1 C 2 C 3 + C 1 C 2 C 3 +^ ]] 

jp [ u^u^^+^u^Uj +2 u^u 2 u[u 3 + 2 u^u^ +u 2 u 3 u[ 2 -^u^u^ +u 2 u 3 u[ 2 + u 3 u 2 uf]] 
+^-J^p|u 2 u 3 u[ 2 + u/c^ +2u 1 u 3 c 1 c 2 + u 2 u 3 c 2 +2u 1 u 2 C|C 3 +2u 1 c i c 2 c 3 +u 2 c 3 c 2 + 



U 3 C 2 C 1 2 ]] 



3Xj L 



u[ 2 c 2 c 3 +2u[u 2 C!C 3 + 2u[u 3 c 1 c 2 + u 2 u 3 c 2 + c 2 c 3 cf jj 

[p [ ^1 U 3 U 2 2 + U 2 2 U 1 U 3 + 2 U l U 2 U 2 U 3 + 2 U 2 U 3 U 1 U 2 + U l U 3 ^2 + 2ll 2 U 1 U 2 U 3 + U l U 3 U 2 2 + U 3 ]] 



^ r~r„/„/„/2 , — t tt 



+- — p u[u 3 u 2 + u 2 CjC 3 +2u 2 u 3 c 1 c 2 +^^02 H^u^c^H^u^c^ +Uj c 3 c 2 + 
3x 2 L L 

+ g|"[p[u 2 2 c^ + 2uXc^+ 2u^c 1 c 2 + u^c 2 +c 1 c 3 c 2 ]J 



^C 2 ]] 



3x. 



p [ u 1 u 2 u 3 2 +u 3 / u[u 2 +2 u 2 u 3 u[u 3 +2 u^u^ +u x u 2 u 3 2 + 2u 3 u[u 2 u 3 +U 2 u i U 3 2 + u i U 2 U 3 2 ]] 

lC 2 + ^c 2 c 2 ]j 



+ - 1 p| U[U 2 U 3 + U 3 0 { C 2 +2U 1 U 3 C 2 C 3 + UjU 2 C 3 +2U2U3CJC3+2U3CJC2C3 + U 2 C 



3x 3 L L 



= -pv 



ufcfz +2u 2 u 3 CjC 3 + 2u|u 3 c 2 c 3 + <u 2 c 2 + c 1 c 2 c 3 |J 
u^c^ +u 3 c 1 c 2 +|c 1 c 2 c 3 j 



^1^2 C 3 + 
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Equation Set 12 



Directional Turbulent Energy (12-1) 




Turbulent Shear (12-2) 
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Equation Set 12 



Directional Turbulent Energy Fluxes (12-3) 
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Equation Set 12 



Directional Turbulent Energy Fluxes (12-4) 



_a_ 

3t 



- a 



+ u, 



+ 2 



UjU^ t^-+u;u 



[«]+<^-|u;ufi+ u 

2 , „'„'2 , ^ U 2 



3 



2"1 



3x t 



+ u 1 u 2 u 3 



+ u 2 u 2 



^2 5u i ^7^772 



■ + u 1 u 1 



+ u : X 2 ^ L 



- 2u;u; 



- u 



+ 2 



/2 



3x 2 3x 
3x 2 L 1 ZJ 9xj I- 1 -I 8x 3 



3xj 



+ Cl c 2 A[ u;2 ] + c ^[ uf ] + CiC3 _^[ u -] 

"77T7 3 



+ 2 



«2 2 ^-[<^]+<u^ ^-[u;u;]+u; u ; ^-[ u ^] 



ax. 



+ <u; — c* + u f — 

dx 2 1 1 ax. 



+ u x u 3 - 



ax, 



u^ + ^u^J + ^u^u? 



uXu--3uX[u^ 2 ]]i^ 



+ 
= 0 



^ 11 ^ _5_P 

p ax 2 

<u]uf - <ul [<] - 2 u]ul [uX ]] A 



u;x 2 -< 2 [uf]-2<ul[u^]]i^ 
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dXj dx 2 dx 3 



dx~ dx 



^72 9u 2 



dx! dx 3 dx 



K U 3] + ^-R U 3] + ^-[ U 



u[u 2 



/2 



dx 



3x 0 L 1 JJ 3x 1 3 



3 

+ u ^ 2 ^[^] +u i u 2^-[c^]+ u y3^-[ c 2c 3 ] 



+ uy 2 — H + u 2 2 — h] + u;u;-[c,c 3 ] 



+ ^[ U 2 U 3] + C 1 C 2 ^[ U 2 U s] + C 1 C 3 ^"N U 3] 



+ CjC 3 



+ Cl c 2 



dXj 

_a_ 

9xj 



+ c 2 c 3 — - 

OX; 



■[<u 2 ] + c^-?-[<u 2 ] 



dx, 



+ c 2 c 3 

dx. 



+ 



+ 



dx 



3x^ 



+ 



u^u; 2 -u^ [u; 2 ]-2<u 2 [u^]]i^- 



Uj'u^Uj 2 - UjU2 



= 0 
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Equation Set 13 (13-1; 13-2; 13-3; 13-4) 
Wi =3c^[cf] 

C ? C 2 =C f [ C 2] + 2C 1 C 2 [ C 1 C 2 ] 

c,c 2 C3 = [cj] + 2c^ [c 2 c 3 ] 
neglected terms being of Burnett order. 

Equation Set 14 ( 14-1; 14-2; 14-3; 14-4) 

uT^[< = 5<ul[uf] 

tftf = + 2 <ul[ u]uT] 

< u 2 uf =uT^"[uf] + 2 ^[ulul] 

neglected terms being of Burnett order. 
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Equation Set 15 

Time Average Continuity (15-1) 
Time Average Momentum (15-2) 

= 0 
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Time Average Directional Thermal Energy (15-3) 



at 



+ 2 
1 

P 

= V 



1 dx x 



+ iu 



dx 0 



c: — — + cx. 



3u, 



3xj 1 2 3x 2 1 5 dx 3 



--hex 



3il 



dx. 



BxJ^ 1 x ] 9x 2 
1= 

-c -cr 
3 *J 









pc 2 cf 


3x 3 


P c 3 cf ] 



Time Average Thermal Shear (15-4) 
|[— ] + 8 [ c =] + ^_L[— ] + 5; 3 [o— ] 



2 

+ c t c 2 ^ + c 2 ^T" + c 2 c 3 — 



+ c 



2 du 2 



3u 1 + ^i 9ui 
3xj 2 3x 

3u 



9xj 



2 = ^ U 2 







3x 2 


P¥2 



+ 



dx. 



\P C i C 2 C 3 ] 



= -vc x c 2 



FILE NO. 004210.P001 



-46- 



GWHVDLM 



Equation Set 15 

Closed Time Average Directional Thermal Energy Fluxes (15-5) 



ar=f 



C,cf 



+ U 2 ~ 

OX 2 L 



+ u. 



dx. 



erf 



+ 3 

+ 3 

- £ 

~ 2 



=f 3uj =f ==f 3u t 
CiCi _ + C2Ci — + c 3Cl — 



dx. 1 





— <2 




— d 








«x 2 


7 


<2x 3 







c t c - 3 CjCj 



(15-6) 
3t 



+ -_a_ 

1 dx. 



CjC 2 



- a 

+ u, 



dx 2 . 



c,c 2 



dx„ 



=f 3u, =f 3uj = j 3uj 

+ C 1 C 2 C 2 C 2 ^ C 3 C 2 ^ 

dXj dx 2 dx 3 



+ 2 



+ 2 



du, 



— 2 5u 2 — j 9u 2 = 

C 2 Cj — h CjC 2 — V CjCjCj — 

dXj dx 2 dx 3 



c,c 2 ^-[c,c 2 ] + c^ ^-[ciC 2 ] + c 2 c 3 ^-[ciC 2 ] 



= |[Ic 2 c*-3c 2 c>] 
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Closed Time Average Thermal Structure (15-7) 



|[ Cl c 2 c 3 ] + Ul A[ CiC2 e 3 ] + A[ CiC2C3 ] + U3 _| 

3x 1 



+ CjC2C 3 + C 3 C 2 + C 2 C 3 V 

dx 2 dx 3 
=r 3u 9 ===== 3u 9 == f du 9 

+ C 3 C 1 ^T + Wa ^ + C 1°3 



3x x 



3x 9 



=j 3u 3 =f 8u 3 8u 3 

+ C 2 C 1 — h ~ H 0^203 — 

dXj dx 2 dx 3 



+ cf ^-[c 2 c 3 ] + c lC2 ^-[c^J + cft ^[c 2 c 3 ] 



+ C 1 C 2 7^[ C 1 C 3] + C 2 ^-[ C 1 C 3] + C 2 C 3 ^"[ C 1 C 3] 



dK 3 



9 [ciCaJ+^Cj^-^J + c^ 



+ CjC 3 



0^! ' 



3 - 



u c { c 2 c 3 
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Time Average Directional Turbulent Energy (15-8) 




Time Average Turbulent Shear (15-9) 
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Closed Time Average Turbulent Directional Energy Fluxes (15-10) 




= 0 
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Equation Set 15 



Closed Time Average Turbulent Directional Energy Fluxes (15-11) 
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Equation Set 15 



Closed Time Average Turbulent Structure (15-12) 



dx n 



, , , Bu, , fl du t , , 2 3u, 



+ u ' u ' ax 



3xj 



3u 2 , „/„/„/ ^ U 2 



-h^U^ 



dx 0 



+ u u 



/„/2 3u 



l w 3 



+ u 2 u 1 



^2 9U 3 , 3^3 . 77^7 ^3 



■ + Ul U 2 



+ u 1 u 2 u 3 



3Xj — 3 —3 



dx. 



3xj 1 



+ c 2 c 3 — 
dx, 



-[u^] + c^[ u ^] 



3x, 



+ c,c 2 + c> g^[utt] + c 2 c 3 ^[u;u 3 ] 



3x, L 1 iJ z 3x 2 L 1 JJ 15 8x 3 



+ <u; ^"[ciC 2 ] + u 2 u^ -^-[c,c 2 ] + uf ^-[c^a] 



3x 0 



3x. 



3x x 1 



0 



FILE NO. 004210.P001 



-52- 



GWHVDLM 



Equation Set 16 



(16-1; 16-2; 16-3) 



c[c™ +1 c n } = /[c^] c|-^ + mp] cjcr^+n^] cjc^T' 
neglected terms being of Burnett order. 



Equation Set 17 

(17-1; 17-2; 17-3) 

< T7 <v=/p] <-^>y+mpiui] «"- 7 <+"[<^] <<X 

V = l pR] u^VV+mp] <uT^ + n[uluI] ufaX 
<<^ = /[uX] <- ; u^u; n +m[ulul] uf^V + np] « m < 

neglected terms being of Burnett order. 
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Equation 18 
(18-1; 18-2; 18-3) 



l+I n m n 
C l C 2 C 3 



/ m+2 n 
C 1 C 2 C 3 



I m n+I 
^lW U 3 



= / 

= 0 

= m 

= 0 

= n 

= 0 



C l C 2 C 3' ^ m rt eVen 



/, m and n all odd 



c l c 2 c 3' m / and n even 
/, m and n all odd 



cjc™ C3 , n odd, / and m even 
/, m and n all odd 



neglected terms being of Burnett order, where /+m+n = 2N +l,an odd number. 
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CLAIMS 



WHAT IS CLAIMED IS: 



1 1 . A method of analyzing anisotropic turbulent flows of an anisotropic fluid 

2 comprising: 

3 defining a set of moment equations governing time average thermal and turbulent 

4 motion, directional kinetic energy, shear, directional kinetic energy fluxes, and structure 

5 correlations; and 

6 defining n m order directional kinetic energy fluxes and structure correlation 

7 equations closure relationships using (n + l) th order density gradient independent time 

8 average thermal and turbulent moment closure relationships to yield a set of closed time 

9 average turbulent moment equations. 



1 2. The method of claim 1 wherein the set of moment equations governing time average 

2 turbulent directional kinetic energy, shear, directional kinetic energy fluxes, and structure 

3 correlations is defined by Equation Set 12. 

1 3. The method of claim 1 wherein n is odd. 

1 4. The method of claim 3 wherein the density gradient independent time average 

2 thermal moment closure relationships are defined by Equation Set 16. 

1 5 . The method of claim 3 wherein the density gradient independent time average 

2 turbulent moment closure relationships are defined by Equation Set 17. 
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1 6. The method of claim 1 further comprising solving (n + 1)* and (n + 2)* order 

2 moment sets. 

1 7. The method! of claim 1 wherein the set of closed time average turbulent moment 

2 equations are defined by Equation Set 15. 

1 8. A method of analyzing time average directional thermal energy in turbulent flows of 

2 an anisotropic fluid by solving Equation 15-3. 

1 9. A method of analyzing time average thermal shear in turbulent flows of an 

2 anisotropic fluid by solving Equation 15-4. 

1 10. A method of analyzing time average directional thermal energy fluxes in turbulent 

2 flows of an anisotropic fluid by solving Equations 15-5 and 15-6. 

1 11. A method of analyzing time average thermal structure correlation in turbulent flows 

2 of an anisotropic fluid by solving Equation 15-7. 

1 12. A method of analyzing time average directional turbulent energy in turbulent flows 

2 of an anisotropic fluid by solving Equation 15-8. 

1 13. A method of analyzing time average turbulent shear in turbulent flows of an 

2 anisotropic fluid by solving Equation 15-9. 

1 14. A method of analyzing time average turbulent directional energy fluxes in turbulent 

2 flows of an anisotropic fluid by solving Equations 15-10 and 15-1 1. 
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1 15. A method of analyzing time average turbulent structure correlation in turbulent flows 

2 of an anisotropic fluid by solving Equation 15-12. 

1 16. A method of analyzing turbulent flows of an isotropic liquid comprising: 

2 defining a set of moment equations governing time average directional kinetic 

3 energy, shear, directional kinetic energy fluxes, and structure correlations; 

4 defining n m order directional kinetic energy fluxes and structure correlation 

5 equations closure relationships using (n + 1)* order density gradient independent time 

6 average thermal and turbulent moment closure relationships to yield a set of closed time 

7 average turbulent moment equations ; 

8 setting all directional thermal energies equal and solving the total thermal energy 

9 equation; 

1 0 setting density equal to a constant. 

1 17. The method of claim 16 further comprising solving the resultant equation set. 

1 18. The method of claim 16 further comprising adding the resulting turbulent flow 

2 equation set to conventional Navier Stokes equations for isotropic fluids and solving the now closed 

3 turbulent Navier Stokes set. 
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ABSTRACT 



A general, closed, anisotropic kinetic turbulence theory for gases and liquids is based on 
new solutions of the Maxwell moment equations of the Boltzmann equations. These solutions 
provide a closed initial equation set for the four time average fluid mechanic variables, the sixteen 
time average thermal motion correlation and the sixteen time average turbulent motion correlations 
listed in Table I. 
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